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Abstract. We analyze the properties of quasar variability using repeated SDSS imaging data 
in five UV-to-far red photometric bands, accurate to 0.02 mag, for ~13,000 spectroscopically 
confirmed quasars. The observed time lags span the range from 3 hours to over 3 years, and 
constrain the quasar variability for rest-frame time lags of up to two years, and at rest-frame 
wavelengths from 1000A to 6000A. We demonstrate that ~66,000 SDSS measurements of mag- 
nitude differences can be described within the measurement noise by a simple function of only 
three free parameters. The addition of POSS data constrains the long-term behavior of quasar 
variability and provides evidence for a turn-over in the structure function. This turn-over indi- 
cates that the characteristic time scale for optical variability of quasars is of the order 1 year. 



1. Rest-frame Time Lag and Wavelength, and Luminosity Dependence 

Significant progress in the description of quasar variability has been recently made by 
employing SDSS data. Vanden Berk et al. (2004, hereafter VB) used a combination of 
imaging and spectrophotometric magnitudes to investigate the correlations of variability 
with rest frame time lag (up to 2 years), luminosity, rest wavelength, redshift, the pres- 
ence of radio and X-ray emission, and the presence of broad absorption line outflows. 
Variability on longer time scales was studied by de Vries, Becker & White (2003, here- 
after dVBW) who compared SDSS and POSS measurements. We extend these studies 
by analyzing repeated SDSS imaging data. Here we summarize the main results obtained 
for a sample of ^13,000 quasars with ~66,000 magnitude difference measurements in 5 
SDSS bandsf (for a detailed description of SDSS survey see Stoughton et al. 2001, and 
references therein), and recalibrated POSS I and POSS II photometry (Sesar et al. 2004). 
Details and supporting analysis will be presented elsewhere (Ivezic et al. 2004, in prep). 

We find that SDSS magnitude difference measurements, Am, follow an exponential 
distribution, p(Am) cx exp(— |Am|/A c ), where the characteristic variability scale, A c 
is a function of rest-frame time lag (At RF , days), wavelength (X RF , A), and absolute 
magnitude in i band (Mi). The variability scale is related to the more commonly used 
structure function (see dVBW for definition) by SF = \/2A c ; we find that the function 



describes Am measurements to within the measurement noise (~0.02 mag). Note that 
there is no dependence on redshift (see Figure 1). The difference between the data and this 
model (right panel) reveals a feature at ^2800 A, which is coincident with the Mgll line. 
This anti-correlation between the continuum and Mg II line variability is indeed known 
from reverberation mapping - due to the excellent SDSS photometry, we recovered this 
result using broad-band data! 

f We do not employ spectrophotometric magnitudes as they are 3-4 times less accurate than 
imaging magnitudes, and are not available for u and z bands. 
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Figure 1. The left panel displays the measured structure function as a function of rest-frame 
time lag and wavelength (each pixel contains ~ 25 objects). The middle panel shows the struc- 
ture function in a narrow range of rest-frame time lag and wavelength (100 d< At^F <150 d, 
1000A< \ RF <2000 A), as a function of redshift and luminosity. The lines of constant vari- 
ability are nearly parallel to the redshift axis, suggesting that the dependence of variability 
on luminosity is much stronger than the dependence on redshift. The right panel displays the 
difference between the data shown in the left panel and the best-fit model described by eq.(l.l). 
The red line shows the composite quasar spectrum from Vanden Berk et al. (2001). The negative 
residuals at ~2800 A are due to the anti-correlation of continuum and Mgll line variability. 
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Figure 2. The dependence of structure function on rest-frame time lag, in the wavelength range 
2000-3000 A, for two data sets: SDSS-SDSS for short time lags (small symbols), and SDSS-POSS 
for long time lags (large symbols). 



2. Detection of a Turn-over in the Structure Function 

The time lags between Palomar Observatory Sky Surveys (POSS I and POSS II) 
and SDSS are much longer (up to ~50 years), than spanned by the available SDSS 
repeated imaging data, and thus offer a possibility of detecting deviations from the simple 
power-law measured for short time scales using repeated SDSS imaging. The observed 
SDSS-POSS long-term variability is smaller than predicted by the extrapolation of the 
power-law from eq. 1.1 (Sesar et al. 2004). The observations (Fig. 2) are well described by 
SF(At RF ) = D (1 - e -( A Wrr) (2.i) 

with D = 0.32 ± 0.03, r = (390 ± 80) days, and 7 = 0.55 ± 0.05. This best-fit is shown 
in Figure 2 by the dot-dashed line. We conclude that the characteristic time scale for 
optical variability of quasars is of the order 1 year. 
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